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ABSTRACT 

The potential of the nonaxisymmetric magnetic instabihty to transport angular momentum and 
to mix chemicals is probed considering the stability of a nearly uniform toroidal field between 
conducting cylinders with different rotation rates. The fluid between the cylinders is assumed 
as incompressible and to be of uniform density. With a linear theory the neutral-stability maps 
for TO = 1 are computed. Rigid rotation must be subAlfvenic to allow instability while for 
differential rotation with negative shear also an unstable domain with superAlfvenic rotation 
exists. The rotational quenching of the magnetic instability is strongest for magnetic Prandtl 
number Pm — 1 and becomes much weaker for Pm 7^ 1. 

The effective angular momentum transport by the instability is directed out- 
wards(inwards) for subrotation(superrotation). The resulting magnetic-induced eddy viscosi- 
ties exceed the microscopic values by factors of 10-100. This is only true for superAlfvenic 
flows; in the strong-field limit the values remain much smaller. 

The same instability also quenches concentration gradients of chemicals by its nonmag- 
netic fluctuations. The corresponding diffusion coefficient remains always smaller than the 
magnetic-generated eddy viscosity. A Schmidt number of order 30 is found as the ratio of 
the effective viscosity and the diffusion coefficient. The magnetic instabihty transports much 
more angular momentum than that it mixes chemicals. 

Key words: instabilities - magnetic fields - angular momentum transport - turbulent diffu- 
sion. 



1 INTRODUCTION 

The current-driven magnetic instability is an important instability 
of to roidal fields w hich is basically no naxisymmetric ('kink-type') 
( iTav ler 1957, 1973; Vandakurovl 1972h . In real fluids a toroidal field 
becomes unstable at a certain magnetic field amplitude depending 
on the radial profile of the field. A global solid- body rotation of 
the system influence s the critical field amplitude jAchesorJl 19781 : 
IPitts & Tavlei|[l985l) . It is strongly increased for fast rotation so 
that we shall speak about a 'rotational quenching' of the Tayler 

instability (TI). 

It has been shown by iKitchatinov & Riidiged 1 I2OO8I) that this 
magnetic instability in rotating spheres with an equatorsymmetric 
magnetic field (peaking at the equator) strongly differs for sub- 
Alfvenic rotation (Q < J7a) and superAlfvenic rotation (O > 
J7a)- Here J7a is the Alfven frequency 



weak fields than for strong fields (Fig.[T). For Oa > the growth 
of the instability is very fast, its growth rate is of the order of the 
rotation rate. The values of the microscopic dissipation parame- 
ters (viscosity, thermal diffusion and magnetic resistivity) have a 
very small influence on the stability behavior of strong fields. The 
growth rate runs linearly with the magnetic field amplitude. 

The situation strongly differs for superAlfvenic flow with 
Oa < ^ where instability only exists for finite values of the heat 
conductivity. It does not exist for both x = and x — 00. The 
latter is obvious as for extremely fast thermal diffusion the density 
fluctuations become zero and the buoyancy term disappears. The 
adiabatic case with x = is also almost withou t dens ity fluctua- 
tions so that no extra instability appears CAchesonlll978h . 

It seems as would the lower limit of the magnetic field for 
neutral instability be fixed by the Roberts number 



Oa = 



(1) 



■n 



(2) 



with Q the rotation rate and R* a characteristic radius. The growth 
rates of the instability are many orders of magnitudes smaller for 
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There is a minimum of the field amplitude sufficient for instabil- 
ity which is well represented in Fig.[I]by the curve for q = 2500. 
Note, however, the smallness of the corresponding growth rates. 
The maximum growth rate in the weak-field domain is ~ 10~*i?, 
i.e. the growth time is more than 1000 rotation times. It is an open 
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Figure 1. The growth rates normalized with the rotation rate for the kink- 
type instability of an equatorsymmetric toroidal field in a sphere, here with 
Pm = 5 • 10"'^, (see Kitchatinov & RUdiger 2008). The two red lines are 
for q = (right) and the solar value q = 2500 (left). The blue line gives 
the limits for both g = and x ^ oo. 



and challanging question whether such slow nonaxisymmetric in- 
stabilities are important for the dynamics of the stellar interior. 

The instability for strong fields with (Qa ^ is much faster 
and it grows linearly with the magnetic field independent of the ro- 
tation rate. The instability does not 'feel' both the thermal diffusion 
and the density stratification. With nonlinear simulations we shall 
answer the question how strong the angular momentum transport 
and the turbulent mixing due to the Tayler instability are. The com- 
plex destabilizing influence of finite thermal diffusion (see Fig. 1 in 
Acheson's paper) is ignored in the calculations what makes a dras- 
tic simplification for the theory. Also for reasons of simplicity the 
computations are done in cylindric Taylor-Couette geometry where 
the differential rotation can be fixed by choice of the rotation rates 
of the inner and the outer cylinder. 

For liquid sodium the Roberts number is small, q = 10""'. 
The neglect of the thermal diffusion is thus fully realistic for MHD 
experiments in the laboratory. It is also realistic in stellar phy sics if 
the sce nario of decaying fossil magnetic fields is considered. IMossI 
( |2003|) finds poloidal fields of order 1000 G surviving in the in- 
ner radiation zones after the Hayashi track. Only after a short time, 
differential rotation produces the toroidal field component which 
tends to dominate the poloidal fields. The resulting Maxwell stress 
suppresses the differential rotation after the Alfven travel time of 
the poloidal field unless the toroidal field becomes unstable. It 
seems to be clear that in order to be effective, the growth rate of 
the instability must be of order of the rotation rate rather than or- 
ders of magnitudes smaller. 

The resulting instabilities form a nonaxisymmetric pattern of 
flow and field which transports both angular momentum and chem- 
icals like lithium. We are interested in the possibility to express 
these fluxes by the gradients of mean quantities. The angular mo- 
mentum flux for rigid rotation in our simplified model proves to be 
very small. It is thus reasonable to express the angular momentum 
transport by an eddy viscosity and also to express the flux of 



chemicals by an eddy diffusivity Dt ■ The amplitudes of the trans- 
port parameters are computed by nonlinear numerical simulations. 



2 ASTROPHYSICAL MOTIVATIONS 

There is a variety of astrophysical applications. Only some of them 
are discussed in the present paper. 



2.1 Rigid rotation in the solar interior 

We know from the results of the helioseismology that the solar 
radiative interior rotates rigidly. This is insofar surprising as the 
microscopic viscosity with v <1Q cm^/s of the hot solar plasma 
yields a diffusion time of more than 10^'' yr so that any initial ro- 
tation law (which certainly existed) should be still existing without 
any other influence. One needs effective viscosities of 10^ cm^/s 
to explain the de cay of an in itial rotation law within the lifetime of 
the Sun (see also lYang & B i 2006). Future results of asteroseismol- 
ogy will show how general this statement is. Also the inclusion of 
the angular momentum transport by a poloidal magnetic field sub- 
ject to differential rotation needs for the explanation of the solar 
spin-down the incr ease of the microscopic visco sity by a few or- 
ders of magnitude jRiidiger & Kitchatinov|[l99^ . We have to test 
whether the angular momentum transport by magnetic instabilities 
of fossil magnetic fields is strong enough to produce the quasi-solid 
inner rotation of the Sun. 



2.2 Diffusion of chemical elements 

The lithium at the surface of cool MS stars decays with a timescale 
of 1 Gyr. It is burned at temperatures in excess of 2.6 x lO" K 
which exists about 40.000 km below the base of the convection 
zone. There must be a diffusion process through this layer to the 
location of the burning temperature. Hence, we have to look for a 
slow process which, however, is only one or two orders of magni- 
tude faster than the molecular diffusion. The molecular diffusion 
at the bottom of the solar convection zone results as Dmicro — 30 
cm^/s teames, Charbormeau & M acGregor '1999") which must be 
increased to about 10^ cm^/s as the consequence of a hypothetical 
instability. Note the smallness of this quantity: the plasma veloci- 
ties for such a diffusion coefficient are by many orders of magni- 
tude smaller than the convective velocities of about 100 m/s at the 
bottom of the convection zone. 

There is thus a good motivation to look for a mechanism 
with a diffusion coefficient Dt > Dmicro but forming a turbulent 
Schmidt number with 

Sc=-gL»i (3) 

(see 'Zahn' '1992')^ The expres sions formulated by 
Eggenberger, Maeder & Mevne 3 j2005l) on the basis of the 
'Tayler-Spruit' dynamo scenario lead for a toroidal magnetic 
field of about 100 Gauss to Sc ~ 100 (where the turbulent 
magnetic diffusivity has been used as the value for Dt). This 
is a reasonable value for the Schmidt number. It is doubtful, 
however, that a toroidal field of 100 Gauss is unstable. For the 
upper part of the solar radiative core we found 600 Gauss as the 
minimum amplitude of a toroidal field to become unstable under 
the presence of the a ctual rotation rate, dens i ty str atification and 
diffusion processes jKitchatinov & Rudigerl l2008h . A dynamo 
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action could not be confirmed so far jZahn. Brun & Mathi j[2007l : 
iGellert. Rudiger & Elstnerll2008h . 

Considering the transport processes in the radiative interior of 
massive (15 A/p) MS stars on the basis of ana lytic expressions by 
ISpruitl j2002h . iMaeder & Mevne^ bOOSlllOOSl) computed viscosi- 
ties up to 10^^ cm'^/s for a rotational velocity of 300 km/s. The 
resulting diffusion time is only a few 100 yr An extreme chemi- 
cal mixing is the immediate consequence. With similar expressions 
lYoon. Langer & Norman (2006) present evolutionary models of ro- 
tating stars with different metalicities. 

Angular momentum transporters could be Maxwell stress by 
both poloidal fields {Br) and toroidal fields (B^) which are con- 
nected by the differential rotation. If the Maxwell stress can be 
expressed by an effective viscosity and a very similar expression 
is used for chemical diffusion processes then the material mixing 
becomes so strong that the resulting evolution models are incom- 
patible with observations. Brott et al. ( 2008) demonstrate that the 
corresponding chemical mixing in massive stars must even be ne- 
glected to understand the observations. 

Note that the material mixing only results by the action of the 
kinetic part of the momentum transport tensor rather than by its 
magnetic part so that the diffusion tensor can be approximated by 



{ui{x,t)u'j{x,t)) 



(4) 



(for details see lRudiger & PiDirill200lh . There is no magnetic influ- 
ence in the diffusion equation except the magnetic suppression of 
the correlation tensor of the fluctuations. 

If the correlation time of the instability scales with its growth 
time then the coefficient of the diffusion of chemicals can simply 
be written as 



7 



(5) 



For strong fields one can estimate the growth rate 7 ~ 
n (see below). Inserting numbers of the solar tachocline 
(Dt ~ 10^ cm^/s, n ~ 10~^ s"^) leads to u'r of or- 
der 0.1 cm/s. It is not easy to produce such a weak turbu- 
lence intensity in radial direction. The idea is that the sta- 
ble density stratification below the convection zone suppresses 
the vertical turbulence co mponent so that a highly anisotropic 
turbulence field evolves 1 Vincent. Michaud & Meneguzzil 1 1 9961 : 
iTogue. Lignieres & VincentI |2006|) . It has been shown, however, 
that due to the stellar rotation even a strictly horizontal flow pat- 
tern obtains radial components which are able to transport chemi- 
cals t hrough the layer below the convection zone dRudiger & PipinI 
l200lh . The result is that the Li depletion becomes anticorrelated 
with the (rapid) stellar rotation or - with other words - the fast ro- 
tators should possess the highest surface concentrations of lithium. 
This consequence of the rotational quench ing of the turbulence ha s 
indeed been found in the Pleiades cluster dXing. Shi & WeilboOTh . 
Such a rotation-lithium correlation, however, could also be the con- 
sequence of the rotational quenching of the TI which is considered 
in the present paper 



3 MAGNETIC TAYLOR-COUETTE FLOW 
3.1 The model 

A Taylor-Couette container is considered confining a toroidal mag- 
netic field with given amplitude at the cylinders rotating with dif- 
ferent rotation rates (see Fig.|2ll. In order to simulate the situation 




Figure 2. Two concentric cylinders with radii i^in and i?,out rotate with 
Oin and I?out- Bg, is the given toroidal magnetic field produced by axial 
currents inside and outside the inner cylinder. 



at the bottom of the convection zone (or even at its top) the gap 
between the cylinders is considered as small. The outer radius is 
Rout and the inner radius is Rin- The fluid confined between the 
cylinders is assumed to be incompressible with uniform density 
and dissipative with kinematic viscosity i' and magnetic diffusiv- 
ity 77. Derived from the conservation law of angular momentum the 
angular velocity of the rotation is 



with 



77^ 



where 



V 



Rin 



b = 



".9 -^in-J'^in 



(6) 



(7) 



(8) 



Rout ^in 

J?in and Qout are the imposed rotation rates of the inner and outer 
cylinders. After the Rayleigh stability criterion the flow is hydrody- 
namically stable for fin > fj^, i.e. /in > 0.25 for fj = 0.5. We are 
only interested in hydrodynamically stable regimes so that fin > 
ff'^ should always be fulfilled. Rotation laws with dO/dR < are 
described by /in < 1; /in = 1 gives rigid rotation. 

Also the magnetic profiles are restricted. The solution of the 
stationary induction equation without inducing flow reads 

B4. = AR+^. (9) 

(in cylindric geometry. A and B are the fundamental quantities; A ■ 
R corresponds to a uniform axial current (7 — 2A/fio) everywhere 
within R < Rout, and B/R corresponds to an additional current 
only within R < Rin - The fields with ^ = are called the current- 
free solutions (current-free between the cylinders). Both terms in 
([9]l are responsible for their own instability where one of them only 
exists in combination with differential rotation. 

In analogy with fin it is useful to define the quantity 

Sout 



flB 



Bin ' 



(10) 



measuring the variation in B^ across the gap. Here we shall con- 
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sider mainly fields which are rather uniform in radius, i.e. /ib ~ 1. 
It is i? = RinRoutA for this choice so that both parts of l|9j are in- 
volved. The condition ^ = leads to profiles with /ib = 0.5, 
and the condition B = leads to ~ 2 both for our stan- 
dard case fi = 0.5. After the stability criterion for toroidal fields 
against axisymmetric perturbations one finds stability of profiles 
with ^5 fiB ^ i/fj (see Shalybkov 2006). The magnetic fields 
discussed in the present paper are thus stable against disturbances 
with m = (see below). 

In the following we shall fix /is = 1 but the rotation law and 
the magnetic Prandtl number 

Pm=-. (11) 
V 

are varied. If the results shall be applied to stellar convection zones 
or e.g. to the fields in galaxies the magnetic Prandtl number must 
be replaced by its turbulent value. 
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Figure 3. The lines of neutral in.stability for A = (AMRI). The curves are 
marked with their magnetic Prandtl numbers. The toroidal fields are current- 
free between the (conducting) cylinders. It is fis = f) = 0.5, fiQ = 0.5. 



3.2 The equations 

The dimensionless incompressible MHD equations are 



Re-^ + Re(u ■ V)u = -VP 



Au + Ha rots x B, 



Rm-^ = AS +Rmrot(u X B), (12) 
with div u — div 5 = and with the Hartmann number 
Ha=4i^. (13) 



Ro = \/Rin{Rout — Rin) is used as the unit of length, rj/Ro as 
the unit of velocity and Bin as the unit of the magnetic fields. Fre- 
quencies including the angular velocity J7 of the rotation are nor- 
malized with the rotation rate Qin- The Reynolds number Re is 
defined as 



Re 



(14) 



and the magnetic Reynolds number as Rm — Pm • Re. The 
Lundquist number is simply S = HaVPm. 

The boundaries are assumed as impenetrable and stress-free 
and as perfect conductors. 

The solutions of the linear equations are free to an arbitrary 
real parameter of any sign. We thus carmot know the sign of the 
flow and/or the field. However, for quadratic expressions such as 
the correlation tensor or the electromotive force we can compute 
the sign as all the solutions are multiplied with one and the same 
parameter. 

Let us apply this idea to the angular momentum transport 



B'rB'^). 

flop 



(15) 



The averaging procedure is thought to be an integration over the 
azimuth (j>. The question we shall answer is whether Tr > for 
dn/dR < and Tfl < for dil/dR > 0. If this is true then 
the angular momentum flows towards the minimum of the angular 
velocity and one can introduce an eddy viscosity i^t, i.e. 

Tr = -j^tR^ (16) 
aR 

with positive fT ■ The sign of vt can be computed with the linear 
theory. We have computed this quantity with both linear and nonlin- 
ear approximations. With the linear theory bifurcation diagrams for 
various rotation laws and magnetic Prandtl numbers are established 



and it is shown that the angular momentum (kinetic plus magnetic) 
flows opposite to the gradient of the angular velocity. It is thus pos- 
sible in this case to define via ( I16t an eddy viscosity on the basis of 
the Tayler instability. With the nonlinear simulations the amplitude 
of this eddy viscosity is computed. 

With U = (ur) the ratio of i/t and Dt results to 



Sc : 



ITr 
q U 



(17) 



with q — — dlogJ?/dlog_R. Without Maxwell stress this quantity is 
certainly of order unity but it can acquire much higher values if the 
magnetic-induced angular momentum transport dominates. 



4 MARGINAL INSTABILITY 

In the linear approximation the solutions can be found in form of 
normal modes in accord to 



F = F{R)e' 



(18) 



Here m is the azimuthal quantum number of the resulting cell pat- 
tern. For instabilities the growth rate 7 — —Q{u>) must be positive 
and the neutral lines are defined by 7 = 0. In this Section the neu- 
tral lines for the magnetic instability with m — 1 are presented. The 
wave number k is optimized so that the resulting Reynolds number 
becomes a minimum for given Hartmann number. 

4.1 The maps of neutral instability 

Let us start with the phenomenon that even current-free toroidal 
fields can become unstable if they are subject to differential rota- 
tion. For = 0.5 follows fiB = 0.5 for the magnetic profile, and 
the differential rotation may be as weak as /in = 0.5 simulating a 
(galactic) rotation law Q oc 1/R. The map for neutral instability 
of this interaction for various magnetic Prandtl numbers is given 
in Fig. [3] The disturbances with m = 1 can grow although the 
magnetic field or the differential rotation alone are stable. We have 
called this phen omenon as Azimuth al Magnetorotational Instabil- 
ity (AMRI, see Rudiger et al."2007, also Ogilvie & Pringle 1996). 
It is very characteristic for AMRI that it disappears for rigid rota- 
tion but also if the (differential) rotation becomes too strong. This 
is the reason why both lines of the marginal stability have the same 
positive slope. 

An axial current (J = 2 A) is not necessary for AMRI. But the 
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Figure 4. The marginal-instability map for toroidal fields with A ~ B 
{fig = 1). Top: Rotation vs magnetic field. The AMRI domain and the 
TI domain are separated by a stable branch. Bottom: The stable branch 
between the AMRI domain and the TI domain coiTesponds to the line 
Mm = 1. It is Pm = 1, /^q = 0.5, fi = 0.5. 



AMRI survives the addition of axial electric currents in the fluid 
conductor (i.e. A =^ 0). It also appears in the instability maps of 
toroidal fields with electric currents between the cylinders (Fig.|4l(. 
In this figure the results for Pm — 1 and for weak subrotation 
(Mn = 0.5) are given for almost homogeneous fields with fiB = 
1. For this magnetic profile both the coefficients A and B in Eq. 
([9} are finite. The characteristic AMRI domain only vanishes for 
magnetic profiles with B — 0. 

Without rotation the magnetic field alone is unstable for Ha ^ 
150 because of the action of the axial current. We shall call this 
domain of subAlfvenic rotation as the domain of the Tayler insta- 
bility (TI). It also exists for B — 0. The critical Hartmann number 
for vanishing rotation does not depend on Pm. In Fig. |4] two un- 
stable domains with different growth rates exist. The growth rate 
is of order Q in the TI domain but it is smaller in the AMRI do- 
main (see below). Both domains are separated by a stability branch 
alon g the line Rm = Ha which disappears, however, for Pm > 1 
(see RUdiger et alj|2007l) . The AMRI domain in Fig. [4] of course, 
disappears for solid-body rotation. 

Note that the rotation basically stabilizes the instability 
jPitts & Tavlei|[l98^ . The faster the rotation the stronger fields 
remain stable. We must keep this effect in mind for stellar appli- 
cations. Magnetic fields which are stable for fast rotation become 
unstable during the stellar evolution with its continuous spin-down. 
Also the magnetic Prandtl number influences the rotational quench- 
ing. One finds the main result already by comparing the instability 
maps for rigid rotation only. The magnetic Prandtl number has been 
varied by two orders of magnitudes in Fig.|5] It makes sense to in- 



Figure 5. Top: The rotational quenching of the TI for rigid rotation and for 
various Pm. The dashed line gives Mm = 1. Note that for Pm = 1 the 
rotational quenching of the magnetic instability is strongest. Bottom: TI is 
a subAlfvenic phenomenon, t} = 0.5, fig = 1- 



terprete the results by means of the 'mixed' Reynolds number 



Rem = \/Re • Rm, 



(19) 



which by definition is symmetric in v and rj as it is the Hartmann 
number. All fluids with the same product v ■ rj have the same Rem 
and Ha. Consequently, the ratio of the mixed Reynolds number 
Rem and the Hartmann number Ha is a magnetic Mach number 

Rem 



Mm = 



Ha ' 



(20) 



in which the numerical values of viscosity and magnetic diffusivity 
do no longer appear. Mm > 1 means superAlfvenic rotation and 
Mm < 1 describes subAlfvenic rotation. 

The rotational stabilization of the toroidal fields also depends 
on the form of the rotation law. It is weaker for subrotation and it is 
stronger for rigid rotation and for superrotation (Fig.|6j. The trend 
is understandable as superrotation is the most stable rotation law 
in hydrodynamics. For given rotation rate the maximum field am- 
plitudes which remain stable are much weaker for subrotation than 
for rigid rotation. With other words: rigid rotation stabilizes mag- 
netic fields more effective than subrotation. For very fast rotators, 
however, the nonaxisymmetric instability is always suppressed. For 
superAlfvenic rotation (Re 3> Ha) the magnetic fields become sta- 
ble also for rotation laws with negative shear. 



4.2 Growth rates 

The growth rates 7 are the negative imaginary part of the frequency 
u! in Eq. l |18t . They are given (normalized with the rotation fre- 
quency of the inner cylinder) in the Fig. [7] The growth time in units 
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Figure 6. The marginal-instability map for subrotation (^n = 0.5), for 
rigid rotation (/in = 1) and for superrotation (/in = 1-25). The curve for 
/ifi = 0.5 is taken from Fig.|4] It is /tg = 1, Pm = 1, f/ = 0.5. 
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Figure 7. The (normalized) growth rates for Pm = 1 (see Fig. |4). Top: 
the TI domain. The curves are marked with their Reynolds numbers. After 
multiplication with Re they are almost identical. Bottom: along the line 
Re = 500 of Fig.|4](AMRI and TI). 

of the rotation time can be obtained as Tgrowth/Trot ~ l/(27r7). 
Hence, for 7 — 1/2tt — 0.16 the growth time equals the rotation 
time. 

The growth rates in the TI domain of Fig. |4] are given in Fig. 
[T] They are much larger than 0.16 so that the growth time of the 
instability is smaller than the rotation period. Note the rotational 
quenching of the growth rates. As they are normalized with the 
inner rotation frequency, however, this is a formal effect. It means 
in Fig. |7] that the physical growth time (in units of the dissipation 
time) does not depend on the rotation rate but it grows linearly with 
the magnetic field (see Goosens, Biront & Tayler 1981). 

Figure|7](bottom) provides the comparison of the growth rates 



of TI with those of AMRI. The values are calculated along the line 
Re = 500 of Fig.m While the growth time of the AMRI remains 
slower than the rotation period the TI modes become faster and 
faster for growing magnetic amplitudes. 



4.3 Magnetic limits 

Figure [5] shows the limit of the subAlfvenic Tayler instability for 
rigid rotation but for various Pm in the Rem-Ha plane. One finds the 
Pm-dependence of the resulting magnetic Mach numbers as rather 
weak. For Pm > 1 the Rem are largest, i.e. for large Pm the ro- 
tational quenching is weak. It is stronger for Pm < 1 but it is 
strongest for Pm — 1 which is the favored application for all kinds 
of numerical simulations. 

The stabilization of the toroidal magnetic fields against Tayler 
instability in rigidly rotating radiative stellar zones is now consid- 
ered. From Fig.|5]we take the simple relation 

Mm ~ 0.5. (21) 

Large Mm means weak rotational quenching of the instability 
which appears to be realized for Pm 7^ 1. If ( 121b is true then 

> 2^/iop u* (22) 

leads to instability, i.e. 

[Gauss] > 2u* [cm/s] (23) 

when the linear rotation velocity u* of the stellar interior is taken 
in cm/s. One thus needs more than 100 kGauss to get a marginal TI 
instability of toroidal fields in the rotating Sun. For Ap stars with 
their higher rotation velocities and higher mean densities even am- 
plitudes of more than 10*' Gauss are necessary. Due to the stellar ro- 
tati on large-scale toroidal fields rema in stable up to this value. Note 
that lHeger. Wooslev & Spruitl ( l2005h work with toroidal field am- 
plitudes of 10* Gauss. For hot MS stars with their extended radia- 
tive zones and a linear rotation rate of (say) 100 km/s the minimum 
magnetic amplitude for the Tayler instability is 10^ Gauss. This re- 
sult only concerns to the fast Tayler instability with growth times 
shorter than the rotation period. By the action of high thermal dif- 
fusion in the radiative stellar cores also weaker toroidal fields can 
become unstable but with much larger growth times (see Fig. [TJ. 
Note that probably almost all rotating objects belong to the fast- 
rotating class with Mm > 1. For the exceptional case of galaxies 
whose inner toroidal fields can be observed, the relation Mm > 1 
is true. 

More calculations with smaller magnetic Prandtl numbers are 
needed to find the final results for the limits of the Tayler instability 
of toroidal fields in rapidly rotating stars. Of course, we can only 
speculate whether the obtained linear scaling also holds for very 
fast rotation. Already small deviations from the linear dependence 
in ( 121b would provide strongly different magnetic field amplitudes. 
On the other hand, the present results from the global theory of 
magnetic instability of toroidal fields in regions of high thermal 
diffusion find the rotational stabilization of TI as very effective. 
Because of the much weaker rotational quenching seen in Fig. |4] 
for the AMRI this sort of magnetic instability seems to be more 
important for stellar physics. We shall see that indeed this instabil- 
ity transports much more angular momentum than the TI. Note that 
the AMRI only exists under the presence of differential rotation. 

Opposite to stars for galaxies we can measure their interior 
magnetic fields so that the magnetic Mach number is well-known. 
It is much larger than unity so that galaxies should be TI stable. 
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However, because of their very high magnetic Prandtl number to- 
gether with the particular galactic rotation law new calculations are 
needed in order to probe their stability. 



5 THE ANGULAR MOMENTUM TRANSPORT 

iGellert & RudigeJ ^20081) have shown that the angular momentum 
transport anticorrelates with Vi?. A formulation as in ( I16t ('dif- 
fusion approximation') is indeed possible. One can show that the 
quantity Tr vanishes for rigid rotation and it is linear in V/?. It 
makes thus sense to define an eddy viscosity with a sofar unknown 
amplitude. In a linear theory a determination of the numerical value 
of the eddy viscosity is not possible. We thus start with nonlinear 
simulations for a container with i) — 0.5, also to minimize the in- 
fluence of the boundary conditions. 

For the same reason the cylinder is regarded with a periodicity 
in the vertical coordinate z set to 6(-Rout — Rin)- We use the MHD 
Fou rier spectral element code describe d by iFournier et al.l j2004h 
and lOellert. Rudiger & Foumiej ( l2007h . In this approach the solu- 
tion is expanded in M Fourier modes in the azimuthal direction. 
The resulting decomposition is a collection of meridional prob- 
lems, each of w hich is solved using a Legendr e spectral element 
method (see e.g. lDeville. Fischer & Munj|2002h . Either M = 8 or 
M = 16 Fourier modes are used, two or three elements in radius 
and twelve or eighteen elements in axial direction, resp. The poly- 
nomial order is varied between N = 8 and A'^ = 16. With a semi- 
implicit approach consisting of second-order backward differenti- 
ation formula and third order Adams-Bashforth for the nonlinear 
forcing terms time stepping is done with second-order accuracy. At 
the inner and outer wall no-slip and perfect conducting boundary 
conditions are applied. 

The Taylor-Couette flow ^ as initial flow and white noise 
for the initial magnetic field are used. Additionally applied is the 
toroidal field with /is — 1 as external field. The differential rota- 
tion is varied between fin — 0.3 and /in = 1.5. At these parameter 
values the field becomes unstable. In the linear regime only the 
Fourier mode m — 1 grows exponentially. Already after roughly 
one rotation, when nonlinear effects become important, also higher 
modes appear. Though, the mode spectrum drops very quickly with 
increasing m. More than 99% of the energy is contained within the 
first three modes. After only a few rotations all excited modes satu- 
rate and a steady-state is reached. The resulting fields are averaged 
in azimuthal direction and the nonaxisymmetric components are 
regarded as fluctuating quantities. 

In accordance to Eq. J15b the angular momentum transport 
is calculated. Figure [8] shows for subrotation (fin — 0.5, quasi- 
galactic rotation law), fiB = 1, Pm — 1 and a Hartmann number 
of Ha — 250 how Tr depends on the magnetic Reynolds num- 
ber. For Tn its values in the center of the gap and averaged over 
z are used. The angular momentum transport scales with Rm"^ up 
to Rm — 600. The last two points do not follow this rule because 
the upper limit for the instability is approached. The differential 
rotation becomes too strong and suppresses the nonaxisymmetric 
magnetic instability equivalent to the results of the linear analysis 
(Fig.|4l(. Maxima as shown in Fig.[8]basically exist due to the limit 
of fast rotation which always exists for the nonaxisymmetric mag- 
netic instability. In order to find the angular momentum transport 
of kink-type instabilities it is necessary to find this maximum and 
to compute the maxima for increasing magnetic fields (see Fig.ll II 
below). The question is how the maxima scale with increasing Ha. 
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Figure 8. Values of Tr (in the center of the gap) for superAlfvenic rotation 
for r) = 0.5, fin = 0.5, fiB = 1, Pm = 1 and for fixed Ha = 250. 
The fitted pai'abola shows Tn depending on Rm'^ . For too high Rm the 
instability is suppressed so that always a maximum exists for fixed Ha at a 
certain rotation rate. 
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Figure 9. Viscosities ut^/u for subAlfvenic rotation and for fin = 0.5. It 

is fiB = 1, Re = 30, ff = 0.5, Pm = 1. 

If they saturate for large Ha the characteristic angular momentum 
transport is known. 

All the nonlinear simulations demonstrate that the angular mo- 
mentum transport Tr scales linearly with dQ/dR so that a char- 
acteristic eddy viscosity can be introduced. In the following, these 
eddy viscosities are calculated in the two instability domains shown 
in Fig. |4] The domain TI is characterized by subAlfvenic rotation 
and strong fields while in the AMRI domain the rotation is super- 
Alfvenic. We shall find that AMRI produces much higher values of 
eddy viscosity rather than TI. 

5.1 Slow rotation 

We start to consider the instability domain of subAlfvenic rotation. 
In the TI area of Fig.|4]for a small Reynolds number (Re = 30) and 
for various Hartmann numbers exceeding Ha = 180 the torque 
Tr and the eddy viscosity are calculated. The eddy viscosity is 
normalized with the molecular viscosity so that the results are fully 
general. The normalized eddy viscosities prove to be very small. 
Note the saturation for increasing Hartmann numbers (Fig. |9j. We 
found the small values of the eddy viscosity as characteristic for the 
considered combination of strong fields and slow rotation. 

5.2 Rapid rotation 

The eddy viscosity is now calculated as a function of the Reynolds 
number of rotation for the same set of parameters used in Fig. [8] 
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Figure 10. The same as in Fig. [8]but for v^jv. Note the existence of the 
maximum value due to the rotational stabilization. 
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Figure 11. The maximum v^/v for various Hartmann numbers, Pm = 1. 
The curves are marked with the Reynolds number where the eddy viscos- 
ity is largest. The normalized eddy viscosity seems to saturate for strong 
magnetic fields. Note that always Mm > 1. 



All these models belong to the domain AMRI in the instability 
map of Fig. |4] The results are plotted in Fig. [TO] For VT/f a 
scaling w ith Rm'^ is found in opposition to the Rm^^ scaling by 
[Maeder & Mevnell (|2005.). Due to the rotational quenching of the 
instability a maximum exists characteristic for the given magnetic 
field. The resulting viscosities ut/v are much higher than for sub- 
Alfvenic rotation but they do not exceed the value of (say) 25. Fig- 
ure [TT] shows the results of calculations of maximum values for 
vt/v for increasing magnetic amplitudes. Also the Reynolds num- 
bers belonging to the maximum viscosity are given. We find a sat- 
uration of the ratio vt/u for strong magnetic fields. The eddy vis- 
cosity exceeds the microscopic viscosity not more than by a factor 
of 100. The same procedure must be applied to the other possible 
magnetic profiles in order to confirm the existence of the magnetic 
saturation suggested by Fig.[TT] The existence of the maximum of 
vt/v for fixed magnetic field as shown in Fig.[TO]seems to be clear. 
The existence of the saturation for stronger magnetic fields must be 
checked by further calculations. 

We have also computed the expression ^ for turbulent dif- 
fusion with growth rates 7 taken from the initial growth phase of 
the instability, which unfortunately is a very simple approximation. 
Therefore, the results do not have the same accuracy as the viscos- 
ity values. Values are approximately one order of magnitude be- 
low the eddy viscosity leading to high Schmidt numbers Sc ~ 30 
(Fig. I12t . We do not find a strong dependence on the magnetic 
field strength. It seems to be true that the kink-type instability for 
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Figure 12. The Schmidt number Sc 
for Ha = 250 and Pm = 1. 
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toroidal fields produces medium viscosity values by Maxwell stress 
and small diffusivity values by Reynolds stress. So far we cannot 
find a remarkable influence of the amplitude of the toroidal field on 
the resulting Schmidt number. 



6 MAIN RESULTS 

The stability of toroidal magnetic fields is considered in cylindric 
geometry. The toroidal fields can be imagined as resulting from the 
interaction of weak poloidal fields with differential rotation. For 
Rm ^ 1 it makes sense to consider the stability for the MHD 
system under neglect of the poloidal field. In the present paper the 
stability/instability of stationary toroidal fields against nonaxisym- 
metric perturbations with m = 1 under the presence of rigid or 
differential rotation is discussed. 

We have shown that under the influence of a rigid rotation the 
marginal instability of toroidal fields is basically suppressed. Note 
that the suppression is strongest for Pm = 1 and becomes weaker 
for any other value of Pm. For differential rotation the situation is 
more complicated. For superrotation {AO /dR > 0) the suppres- 
sion is even stronger than for rigid rotation, while for subrotation 
(dn/dR < 0) it also exists but only for very fast rotation. These 
results are understandable because too strong differential rotation 
always disturbs nonaxisymmetric magnetic configurations. 

For radiative zones the kink-type instability can be consid- 
ered as an angular momentum transporter. The angular momentum 
transport vanishes for solid-body rotation and the angular momen- 
tum is always transported opposite to the gradient of the angular 
velocity so that a 'diffusion approximation with an eddy diffusiv- 
ity makes sense. The eddy viscosity resulting from our simulations 
exceeds the microscopic viscosity not more than by 1-2 orders of 
magnitudes. The truth of this statement bases on the existence of 
the magnetic saturation of the maximum viscosity values presented 
in Fig.[TT] Further simulations are needed to confirm this finding. 

If angular momentum is transported and the transport is not 
completely magnetic then also the considered fluid is continously 
mixed. For stron g mixing the ste llar evolution would too strongly 
be affected (see Brott et al] 120081 . hence the mixing induced by 
magnetic instability must remain weak. Also the lithium observa- 
tions of solar-type stars are not understandable without the exis- 
tence of a rather mild turbulence beneath the stellar convection 
zones. One may argue that the mixing is due to the Reynolds 
stress while the angular momentum transport is mainly due to 
the Maxwell stress. If the Maxwell stress dominates the Re ynolds 
stress then the general transport problem formulated by IZahnI 
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( Il992h is solved. We indeed find such a tendency with nonlinear 
simulations. The eddy viscosity exceeds the diffusion coefficient 
by one or t wo orders of magnitude. Very similar results are reported 
bvlCarballi do. Stone & Pringle (2005) and Johansen, Kla hr & Med 
( l2006h for simulations of the magnetorotational instability. 

It is an open question whether the dominance of eddy viscos- 
ity over turbulent diffusion can even be stronger. Figure [72l does 
not lead in this direction. More simulations with different magnetic 
profiles and/or with smaller magnetic Prandtl numbers must show 
how many orders the magnetic-induced angular momentum trans- 
port exceeds the diffusion of chemicals. 

Note that the values derived with our model ~ in particular the 
value for the diffusion coefficient l|5j - are maximum values as the 
suppressing action of the density stratification and the destabilizing 
action of the microscopic thermal diffusion is still neglected. Also 
the restriction to magnetic Prandtl number of order unity and the 
consideration of a special magnetic profile are still serious limita- 
tions. 
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